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Figure  Captions 


1.  Excitation  energy  versus  the  dimensionless  parameter  qr^  in  units  of  the 
Rydberg  at  several  values  of  . 

2.  Structure  factor  versus  kr  at  various  values  of  R  . 

o  s 

1.  Structure  factor  versus  kr^ .  The  solid  and  dashed  lines  correspond  to  two 

and  three  dimensions  for  R  -  1  and  3. 

s 

4.  Pa  ir  -  cor  re  lat  ion  function  g(r)  versus  r  for  various  values  of  R  .  The 

s 

solid  and  dashed  lines  correspond  to  two  ar.d  three  dimensions. 


iO 


Table  1.  Ground-scate  energies  obtained  from  various  methods  for  the  two-  and 
three-dimensional  charged-boson  system. 


2D  3D 


R 

s 

-Eq :  SCFA 

.  17 

RDA 

SCFA8 

RDA19 

VCA5 

1 

1.2617 

1.2918 

0.7712 

0.8030 

0.7767 

2 

0.7808 

0.8138 

0.4472 

0.4475 

0.4516 

3 

0.5965 

0.6210 

0.3231 

0.3523 

0.3270 

5 

0.4245 

0.4418 

0.2129 

0.2402 

0.2159 

10 

0.2670 

0.2783 

0.1188 

0.1428 

0.1209 

Table  2. 

Pressure  of  the 

2D  charged-boson 

system  obtained 

from  the 

SCFA  and 

RDA. 


2D 


R 

s 

-P/p:  SCFA 

RDA 

i 

0.4060 

0.4306 

2 

0.2557 

0.2713 

3 

0.1961 

0.2070 

5 

0.1404 

0.1473 

10 


0.0911 


0.0928 
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and  for  the  3D  charged-boson  system  the  ground-state  energy  per  particle  is 
evaluated  by  various  methods : 


EQ  -  -  0.8030  Rs3/4  ,  (RDA) 19  (4.3) 
Eq  -  -  0.8030  Rg3/4  +  0.027  ,  (VCA)5  (4.4) 
Eq  -  -  0.8030  R'3/4  +  0.032  .  (SCFA)8  (4.5) 


The  extra  numerical  terms  in  Eqs .  (4.2),  (4.4)  and  (4.5)  are  due  to  the  better 
estimation  of  the  correlation.  The  numerical  data  for  the  ground-state 
energies  in  the  various  calculations  are  summarized  in  Table  1,  where  we  see 
that  the  SCFA  results  are  an  improvement  over  the  RDA  results. 

Differentiating  Eq.  (3.20)  with  respect  to  Rg ,  we  obtain  the  pressure  of 
the  2D  charged-boson  system: 


P 

P 


R  dEn 
s  _ 0 

2  dR 

s 


(4.6) 


The  numerical  results  for  the  pressure  are  listed  in  Table  2  in  comparison 
with  the  RDA  results. 

In  conclusion,  we  remark  that  the  behaviors  of  the  elementary  excitation 
spectrum,  structure  factor  and  pair-correlation  functions  obtained  from  the 
SCFA,  which  includes  the  short  range  correlations  for  the  2D  charged-boson 
system,  are  very  much  like  those  of  the  3D  case.  The  SCFA  results  for  the 
ground-state  energy  are  an  improvement  over  the  RDA  results. 
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factor  increases  more  slowly  than  the  3D  structure  factor  in  the  short- 
wavelength  region  but  more  quickly  in  the  long-wavelength  region  as  qr  is 
varied. 

With  the  use  of  the  above  calculations  for  the  structure  factor,  we 

evaluate  the  pair- correlation  function  from  Eq.  (3.13).  The  result  is  given 

and  compared  with  the  3D  case  in  Fig.  4.  The  pair-correlation  function  starts 

oscillating  at  large  distances,  which  is  not  displayed  in  the  figures.  These 

oscillations  have  very  small  and  broad  amplitudes.  This  long-distance 

18 

behavior  is  similar  to  the  3D  case.  For  Bose  fluids  interacting  via  a  soft 

potential  with  a  Lennard- Jones  type  tail  and  pseudopotential,  in  the  RDA  the 

pair-distribution  functions  become  negative  at  short  distances  and  decrease  as 
-  3 

r  at  large  distances.  This  decrease  corresponds  to  the  existence  of  a 
phonon  spectrum  for  small  momenta.  In  the  RDA  the  pair-distribution  function 
for  the  2D  charged  boson  becomes  negative  with  decreasing  Rg  and  diverges  at  r 
-*  0.  Consideration  of  the  short-range  correlation  of  the  charged  boson 
through  the  local-field  correction  SCFA  improves  the  result  from  the  RDA, 
i.e.  ,  the  value  of  g(0)  obtained  in  SCFA  is  negative,  but  so  small  that  for 
the  practical  purposes  one  can  consider  g(0)  to  be  zero. 

We  have  evaluated  the  ground- state  energy  by  a  numerical  self-consistent 
solution  through  Eqs.  (3.14)  and  (3.19)  to  give 

E_  -  -  1.2918  R'2/3  +  0.03  .  (4.1) 

u  s 

In  the  RDA  we  have  obtained  the  ground-state  energy  per  particle  in  terms  of 

R  as 
s 


E.  -  -  1.2918  R*2/3 
0  s 


(4.2) 


10 


EQ  -  21/3R‘2/3  j  da  j  dq  [S(a,q)  -  1]  (3.20) 

We  note  that  the  3D  ground-state  energy  per  particle  in  Rydbergs  is  given  by 

E_  -  1  31/4  R‘3/4  f  da  P  dq  [S(a,q)  -  1]  (3.21) 

0  *  3  J0  J0 

4.  Results  and  discussion 

In  the  previous  sections  we  have  evaluated  the  elementary  excitation 

spectrum,  pair - corre lation  functions  and  the  ground-state  energy  of  the 

charged-boson  system  from  the  determination  of  the  structure  factor  in  the 

self-consistent  field  approximation.  Figure  1  illustrates  the  elementary 

excitation  spectrum  as  a  function  of  qr  for  various  values  of  R  .  We  see 

o  s 

that  the  third  term  in  the  bracket  of  Eq .  (3.10)  is  dominant  in  the  high- 
momentum  region,  and  thus  the  excitation  energy  is  almost  identical  with  that 
of  a  free-particle .  As  Rg  increases  the  excitation  energy  reduces  more 
rapidly  to  the  free  particle  case.  We  notice  that  in  the  low-momentum  region 
the  3D  excitation  energy  obtained  from  the  SCFA  decreases  as  q  increases  ,  and 
this  reduction  is  quite  significant  with  increasing  Rg .  This  reduction  does 
not  appear  in  the  2D  case  for  the  range  of  densities  1  <  Rg  <  10.  A 
comparison  of  Eqs.  (3.6)  and  (3.12)  shows  that  the  result  in  the  SCFA  gives  a 
correction  to  Eq.  (3.12)  derived  from  the  RDA. 

In  Fig.  2  we  have  given  the  numerical  calculation  of  Eq.  (3.14)  as  a 
function  of  qr  for  various  values  of  R  by  the  iteration  method.  The 
structure  factor  converges  very  rapidly  to  unity  as  Rg  increases.  We  have 
compared  the  2D  structure  factor  to  the  3D  case  in  Fig.  3.  The  2D  structure 


9 


by  the  method  of  iteration.  The  3D  structure  factor  in  the  long-wavelength 
approximation  is  given  by 


S(q)  -  (1 


12 

Rjq4 


(1 


7q2)  ] 


(3.16) 


with 


3 

2R 

7  -  -  -gj  J  dq[S(q)  -  1]  .  (3.17) 

Comparing  Eqs .  (3.9)  and  Eq.  (3.17),  we  find  that  both  forms  are  similar  to 
each  other,  with  dependence  of  the  spatial  dimensionality  on  the  dimensionless 
parameter  Rg . 

From  the  pair-correlation  function  or  structure  factor,  we  can  calculate 
the  ground-state  energy  of  the  charged-boson  system.  We  may  write  the 
interaction  energy  as 


dr  0(r,a)[g(r)  -  l]r  ,  (3.18) 

or 

Eint(a)  "  2  dq  atS<Q’^  -  11  ■  (319) 

2  k 

Expressing  the  wave  number  and  the  density  in  units  of  (2npaQ)  and  Rg ,  we  can 
write  the  ground-state  energy  per  particle  in  Rydbergs  as 


3 


E(q)  -  JffJd  + 


*2  4 

2  2' 
2a  Q 


(3.12) 


The  short-range  correlation  functions  occurring  in  the  Coulomb 
interaction  between  the  charged  bosons  are  expressed  by  the  pair- correlation 
function  g(r) ,  which  represents  the  probability  of  finding  two  boson  particles 
separated  by  a  distance  r.  The  inverse  Fourier  transform  of  the  structure 
factor  yields 


K  f  - 

g(r)  -  1  +  J  dq  q  JQ(qr)(S(q)  -  1' 


(3.13) 


and  from  Eqs .  (2.4),  (2.7)  and  (3.6)  the  structure  factor  becomes 


S(q)  “  (1  +  “All  *  G(q)n'4 

Rsq 


(3.14) 


with 


R2  rrq 


G(q)  -  J4  i|  dk  lcE(|j)[S(k)  -  1]  +  q  j*  dk( [1  -  (^)2]K(^) 


-hi: 


+  (jj)2  E(]j)}[S(k)  -  1]J>  , 


(3.15) 


where  r,  q  and  k  are  expressed  in  units  of  the  Bohr  radius  aQ  and  aQ  , 
respectively,  and  K(x)  and  E(x)  are  the  complete  elliptic  integrals  of  the 
first  and  second  kinds.  The  numerical  solution  of  Eq.  (3.14)  can  be  obtained 


S (q  -  k)  -  S (q)  -  qcostf  +  ... 


(3.7) 


Che  local  field  correction  G(q)  becomes 

G(q)  -  7q 


(3.8) 


and 


7 


d2q  [ S (q) 


1] 


(3.9) 


Then  the  excitation  spectrum  can  be  expressed  as 


E(q)  -  Hfi[L 


2 

2 


4 

q  + 

8m  n 


where  Q  -  (2rcpe  q/m)  is  the  two-dimensional  plasma  frequency, 
three -dimensions  the  excitation  spectrum  in  the  long-wavelength 
is  given  by 


(3.10) 


In  the  case  of 
approximation 


W2  4 

E(q)  -  [ l  -  J  q  +  22] 

*  2m  w 

P 


(3.11) 


where  -  (4ir pe  /m)  is  the  3D  plasma  frequency.  We  remark  that  one  of  the 
authors  has  obtained  the  excitation  spectrum  of  the  2D  charged-boson  system 
through  the  RDA  in  the  long-wave length  approximation  as^2 


6 


where  A  -  Sxp/a^.  The  total  induced  charge  Q  is 


I 


Q  -  -e  d  r  Sp{ r) 


I 


ze  I  d  q 


-sin  Car,) 


G(q)  -  (J)3  -  1 


-  -  ze 


(3.4) 


Equation  (3.4)  indicates  that  the  charged  impurity  is  completely  screened  at 
long  distances.  However,  the  induced  charge  density  [Eq.  (3.3)]  diverges  at 
r  -  0.  This  divergence  is  due  to  the  fact  that  the  linearized  equation  of 
motion  for  the  classical  one-particle  distribution  function  is  invalid  near 
the  charged  impurity.  The  divergence  can  be  avoided  by  taking  quantum  effects 
into  consideration. 

The  elementary  excitation  spectrum  E(q)  is  determined  from  the  pole  of 
the  density -density  response  function  x(q,t*>),  which  yields 

[w(q)  +  i»7  ] 2  *  «2(q)  -  2pe(q)V>(q)  -  0  ,  (3.5) 


and  thus  the  excitation  energy  E(q)  -  Hw(q)  can  also  be  written  as 


E(q)  -  [«(q)2  +  2pe  (q)^(q)  ]  ** 


(3.6) 


We  notice  that  the  elementary  excitation  spectrum  in  the  ring-diagram 

16  -♦  -*• 
approximation  (RDA)  can  be  obtained  under  the  condition  V’(q)  ”  <£(q)  ,  i.e., 

the  neglection  of  the  local  corrections  in  Eq .  (3.6).  Making  use  of  the 

following  expression  for  the  structure  in  the  long-wavelength  approximation, 


We  find  the  ground-state  energy  to  be 


D 


2 


where  E.^t(a)  Is  the  interaction  energy  as  a  function  of  the  coupling  constant 
a  which  is  a  measure  of  the  strength  of  the  coupling  between  bosons. 


3.  Excitation  Spectrum,  Structure  Factor,  Correlation  Function  and  Ground- 
State  Energy 

We  first  investigate  the  response  of  the  charged-boson  system  to  a 
static  impurity  with  charge  ze  located  at  the  origin,  where  the  external 
potential  is 

,  2 

^ext(q.«)  -  *(*)  ■  (3.D 

The  induced  charged  density,  which  characterizes  the  linear  response  to  an 
external  potential  from  a  fixed  charge,  can  be  written  as 


6p(q,w)  -  -  x(q-w)e^ext(q,w) 


(3.2) 


Through  the  inverse  Fourier  transform  of  Eq .  (3.2),  we  obtain  the  induced 

— ► 

charge  density  at  position  r  as 


Sp(  r) 


dq 


sin(ar) 


G(q)  - 


'?>3 


-  1 


(3.3) 


**♦ 

In  Eq .  (2.2),  *o(q,a>)  is  che  dens  i  ty  -  de  ns  i  ty  response  function  for  a 

noninteracting  charged-boson  system  at  T  -  0  given  as 


xQ(q.w) 


=  2pe  (q)/[  (<j  +  ip)' 


« (q)2] 


(2.3) 


where  «(q)  -  q  /2m  is  the  free  particle  energy,  p  is  a  positive 

— ¥ 

infinitesimal  quantity,  and  V(q)  is  the  self-consistent  effective  potential, 


V>(q)  -  4>(q)  ti  -  G(q)  ] 


(2.4) 


Here,  <£(q)  -  2ne  /q  is  the  two-dimensional  Fourier  transform  of  the  Coulomb 
2 

interaction  e  /r,  and  G(q)  is  given  by 


g(^}  “  ■  H  [s(q  -  v  -  ^  ■  (2-5> 

P  J  qK  (2 n)1 

In  Eq.  (2.5)  the  static  structure  factor  S(q)  ,  which  is  the  Fourier  transform 
of  the  pair-correlation  function  g(r) ,  can  be  expressed  as 


S(k)  -  1  +  p 


[g(r) 


l]e 


-  ik*  r 


(2.6) 


The  singularities  of  the  density -density  response  function  represent  the 

energies  of  the  excited  states,  and  the  excitation  energy  of  the  system  is 
— ♦ 

related  to  S(q)  through  the  Feynman  expression 


E(q)  -  <(q)/S(q) 


(2.7) 


spectrum,  structure  factor,  pair - corre la t ion  function  and  the  GSE  of  2D 
char  ged  -  b  o s on  systems.  Therefore,  in  this  paper  we  evaluate  the  above 
quantities  of  a  2D  charged  boson  system,  which  consists  of  N  identical  bosons 
with  charge  e  and  mass  m,  interacting  via  a  Coulomb  potential  at  T  -  0  over 
the  range  of  densities  1  <  <  10 .  To  calculate  the  above  quantities  we 

adopt  the  SCFA  given  by  Singwi  et  al.  We  survey  the  basic  formula  in  Sec.  2, 
and  starting  with  the  numerical  self-consistent  evaluation  of  the  static 
structure  factor  S(q),  we  obtain  the  elementary  excitation  spectrum, 
correlation  function  and  ground-state  energy  in  Sec.  3.  Finally,  in  Sec.  U  we 
present  our  numerical  results  for  the  above  quantities  in  comparison  with 
other  works  in  terms  of  graphs  and  tables. 

2.  Basic  Formulas 

The  self-consistent  field  approximation  in  the  formalism  of  Singwi  et  al 
includes  the  decoupling  of  the  two-particle  distribution  function  in  the 
Liouville  equation  into  the  product  of  two  one-particle  distribution  functions 
and  a  pair-correlation  function, 

f(r,p,r' ,p' , t)  -  f1(r,p|t)  f1(r',p'|t)  g(r  -  r')  ,  (2.1) 

-4  —4 

where  r  and  p  are  the  position  and  momentum  of  each  particle  and  g(r) 

represents  the  equilibrium  static  pair-correlation  function.  The  density- 

—¥ 

density  response  function  x(q.w)  in  Fourier  space  for  an  interacting  system 
becomes 

x(q,w)  -  x0(q.w)/(i  -  ’Kq)x0(q.w) ] 


(2.2) 


1.  Introduction 


Since  Foldy's  pioneering  work  on  the  charged-boson  system,  there  has 

been  continuing  interest  in  this  system  from  the  view  point  of  real  physical 

2 

systems.  It  should  be  pointed  out  that  the  special  characteristic  shared  by 

the  charged  boson  and  electron  gas  is  that  the  properties  of  the  ground  state 

can  be  described  by  only  a  single  dimensionless  parameter  R  -  r  /a  where 

s  o  o 

-H 

r  -  (jt p)  is  the  mean  particle  distance,  p  is  the  number  density,  and  a  - 
2  2 

Y  /me  is  the  Bohr  radius  of  the  particle  with  mass  m.  The  electron  gas  has 

been  widely  investigated  for  its  applications  to  metals.  H-.vever,  the 

charged-boson  system  has  been  largely  ignored  because  of  its  nonexistence. 

Concerning  the  ground-state  energy  (GSE)  ,  Foldy  first  calculated  the 

GSE  and  elementary  excitation  spectrum  of  the-  charged  boson  system  in  the  high 

3 

density  region  (Rg  <  1),  which  was  also  later  evaluated  by  othors.  In  the 

4  5  6 

intermediate  density  region  (1  <  Rs  <  100),  Lee,  Lee  and  Ree  and  Monnier 
evaluated  the  GSE  through  the  variational  method  with  a  Jastrow  trial  wave 
function.  More  recently  Hansen  and  Maxighi^  obtained  a  variational  upper 
bound  to  the  GSE  over  a  wide  range  of  densities  (1  <  R^  <  200)  with  the  use  of 
a  variational  Jastrow  trial  wave  function,  the  hypernetted  chain  integral 

g 

equation,  and  the  Monte  Carlo  method,  and  Hipolito  et  al  investigated  the 

dielectric  properties  and  also  the  two-dimensional  (2D)  and  three-dimensional 

(3D)  classical  electron  systems  by  adopting  the  self-consistent  field 

9 

approximation  (SCFA)  introduced  by  Singwi  et  al. 

Although,  in  the  past  two  decades,  significant  progress  has  been  made  in 
the  study  of  the  2D  electron  systems  for  the  dielectric  function,  structure 

factor  and  pair-correlation  function,  ^  GSE  and  specific  heat,^  effective 

12  13  1^  1 5 

mass,  superlattice,  quantized  Hall  effect  and  other  quantities,  there 

is  much  less  information  about  the  properties  of  the  elementary  excitation 
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field  approximation;  R  -  (a  irp)  ,  where  a  is  the  Bohr  radius  and  p  is  the 

so  o 

number  density.  Starting  with  numerical  self-consistent  calculations  of  the 

static  structure  factor,  the  elementary  excitation,  pair-correlation 

functions,  pressure  and  ground-state  energy  are  evaluated.  These  results  are 

compared  with  those  of  the  two-dimensional  and  three-dimensional  systems 

obtained  from  other  methods.  The  ground-state  energy  is  given  as  - 
-  2/3 

-1.2918  R^  +  0.03,  which  improves  the  result  from  the  ring-diagram 
approximation . 
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